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<N ; ABSTRACT 

3 : 

qq ■ We study prompt optical emission from reverse shocks in the wind-type 

gamma-ray bursts. The emission is evaluated in both the thick and thin shell 
regimes. We discuss the angular time delay effect and the post-shock evolution of 
the fireball ejecta, which determine the decay index of the prompt optical emis- 
sion and the duration of the radio flare. We discuss distinct emission signatures 
of the wind environment compared with the constant interstellar medium envi- 
ronment. We also present two recipes for directly constraining the initial Lorentz 
factor of the fireball using the reverse and forward shock optical afterglow data 
o for the wind case. 
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1. Introduction 

The gamma-ray burst (GRB) afterglow observations are usually explained by a model 
in which a relativistic fireball shell (ejecta) is expanding into a uniform interstellar medium 
(ISM) . However, there is observational evidence suggesting a link between GRBs and massive 
stars or star formation (e.g. Meszaros 2001 for a review). An important consequence of a 
massive star origin for the afterglow is that the fireball shell should be expanding into a 
pre-burst stellar wind of the progenitor star with a density distribution of p oc R~ 2 (e.g. 
Chevalier & Li 1999; Meszaros, Rees & Wijers 1998; Dai & Lu 1998). This wind model is 
discussed to be consistent with some GRB afterglow data (Chevalier & Li 1999, 2000). 



It is expected that many early optical afterglows will be discovered soon in the obser- 
vational campaigns led by HETE-2 and Swift. Long time span (from right after the GRB 
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trigger to a year) observational data will allow us to distinguish differences between the ISM 
and wind models clearly. 

In this paper, we will discuss the optical reverse shock emission for the wind model in 
detail. The previous study (Chevalier & Li 2000) gave the discussions for the thick shell 
case, and we will consider both the thin and thick shell cases. We show that the angular 
time delay effect plays an important role in discussing the light curve decaying phase of the 
reverse shock emission. 



2. The Model 

We consider a relativistic shell (fireball ejecta) with an isotropic energy E, an initial 
Lorentz factor rj and an initial width A expanding into a surrounding medium with den- 
sity distribution of p = AR~ 2 (wind model). In this paper, Lorentz factors 7, radii R and 
widths A are measured in the laboratory frame where the progenitor is at rest. Thermody- 
namic quantities: mass densities p and internal energy e are measured in the fluid comoving 
frame. The interaction between the shell and the wind is described by two shocks: a for- 
ward shock propagating into the wind and a reverse shock propagating into the shell. The 
shocks accelerate electrons in the shell and wind material, and the electrons emit photons via 
synchrotron-cyclotron process. For the synchrotron process, the spectrum of each shock emis- 
sion is described by a broken power law with a peak F v max oc N e B^ and break frequencies: a 
typical frequency v m oc -877^ and a cooling frequency v c oc 1/B 3r yt 2 (Sari, Piran & Narayan 
1998) where N e is the number of electrons accelerated by the shock, B is the magnetic field 
strength behind the shock, 7 and 7 m are the bulk Lorentz factor of the shocked material 
and the random Lorentz factor of the typical electrons in the shocked material, respectively. 
Assuming that constant fractions (e e and e#) of the internal energy e produced by the shock 
go into the electrons and the magnetic field, we get v m oc 7p~ 2 e 5//2 , v c oc 7 _1 e~ 3 / 2 t~ 2 and 



3. Forward Shock 

Observations of optical afterglows usually start around several hours after the burst 
trigger. At such a late time, the shocked wind material forms a relativistic blast wave and 
carries almost all the energy of the system. Chevalier & Li (1999) gave the characteristics 
of the forward shock (blast wave) emission as follows, 

u mJ (t) ~ 1.6 x 10 u C 1/2 <-i^2^V /2 Hz, (1) 
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v c , f (t) ~ 6.3 x l^-^\El?A:HT Hz, (2) 
~ 4.5 rf- 2 C 3/2 eB /2 -2^52 /2 ^^ 1/2 mJy, (3) 

where £ = (1 + z)/2, d = cIl(z)/(2 x 10 28 cm), z and cIl(z) are the redshift and luminosity 
distance of the burst, respectively, 6^(1) ~ 2 x 10 28 cm for the standard cosmological pa- 
rameters {Q m = 0.3, Q A = 0.7 and h = 0.7), e e _i = e e /0.1, e Bi _ 2 = e B /0.01, £ 52 = £/10 52 
ergs, A* = A/5 x 10 11 g cm -1 , td is the observer's time in units of days. The optical 
flux from the forward shock decays proportional to t -1 / 4 initially, and decays faster as 
£-(3p-2)/4 a ft er a transition of u m j through the optical band v R ~ 5 x 10 14 Hz (Chevalier 
& Li 2000) where p is the index of the power law distribution of the accelerated electrons. 
Using z/r* = vr/5 x 10 14 Hz, the break time t m j (the passage of v m j) and the optical flux 
(in the fast-cooling regime) at that time are 

f ^ qfl ^V3 4/ 3 „V3 pl/3 -2/3 . 

*m,/ ~ 30 Q e e - 1 e B _ 2 E 52 v Rif mm, (4) 
4. Reverse Shock 

At earlier time when the reverse shock crosses the shell, the forward-shocked wind and 
the reverse shocked shell carry comparable amount of energy. A significant emission is 
expected from the reverse shock also (Meszaros & Rees 1997; Sari & Piran 1999a). During 
the reverse shock crossing, there are four regions separated by the two shocks: the wind 
(denoted by the subscript 1), the shocked wind (2), the shocked shell material (3) and the 
unshocked shell material (4). Using the jump conditions for the shocks and the equality of 
pressure and velocity along the contact discontinuity, we can estimate the Lorentz factor 7j, 
the internal energy and the mass density pi in the shocked regions as functions of three 
variables 74 (= 77), p\ and p 4 (Sari and Piran 1995). 

There are two limits to get a simple analytic solution to the hydrodynamic quantities at 
a shell radius R (Sari and Piran 1995). If the Lorentz factor is low rf <C / where / = Pa/ Pi, 
the reverse shock is Newtonian which means that the Lorentz factor of the shocked shell 
material is almost unity in the frame of the unshocked shell material. It is too weak to slow 
down the shell effectively so that 73 ~ 77. On the other hand, if the Lorentz factor is high 
rf ^> /, the reverse shock is relativistic, and considerably decelerates the shell material, 
hence 73 ~ r/ 1 / 2 / 1 / 4 . 
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4.1. Critical Radii 

Since the density ratio / is generally a function of R, there is a possibility that the 
reverse shock evolves from Newtonian to relativistic during the propagation. The evolution 
of the reverse shock depends on the ratio of two radii: R y = E /4ttAc 2 t) 2 where the forward 
shock sweeps a mass of E/c 2 rf and R s = A r] 2 where the shell begins to spread if the 
initial Lorentz factor varies by order 77 (Sari & Piran 1995; Kobayashi, Piran & Sari 1999). 
Another important radius is R x where the reverse shock crosses the shell. The lab-frame 
time it takes for the reverse shock to cross a width dx of the shell material is given by 
dtiab ~ dR/c ~ rjf l l 2 dx/c (Sari & Piran 1995). We can regard R/c as time t[ a b in the 
laboratory frame because of the highly relativistic expansion of the shell. Since the whole 
shell width is A ~ max[A , R/r] 2 ], we obtain R x ~ max[(i? s i? 7 ) 1 / 2 , i? 7 ]. 

Considering if j j = max[i? s /i? 7 , R/R^], we can classify the evolution of reverse shocks 
into two cases by using a critical Lorentz factor r\ c = (E / 'An Ac 2 Aq) 1 ' ' 4 (see Sari & Piran 1995 
and Kobayashi & Zhang 2003 for the ISM model). If R s /R^ = (r?/^) 4 > 1 (Thick Shell 
Case), the reverse shock is relativistic from the beginning (at the end of the GRB phase), 
which is different from the ISM model in which the reverse shock only becomes relativistic 
later. The reverse shock crosses the shell at R x ~ (i? s i? 7 ) 1//2 before the shell begins to spread 
at R s: and it significantly decelerates the shell material 73 ~ rj c . If R s /R 1 < 1 (Thin Shell 
Case), the reverse shock is initially Newtonian and becomes only mildly relativistic when 
it traverses the shell at R x ~ i? 7 . We can regard 73 as constant (~ rj) during the shock 
crossing. 



4.2. Synchrotron Emission 

Since the wind density at the initial interaction is much larger than the medium density 
for the ISM model, the cooling frequency v c r of the reverse shock emission is much lower 
than the typical (injection) frequency z/ m r in the wind model. The random Lorentz factor 
7 C of the electrons that radiate at the cooling frequency u c , r could be sub-relativistic or 
Newtonian, 7c (t x ) ~ IC^^^A^E^t 3 ^ for our typical parameters. This makes the 
radiation mechanism cyclotron radiation at low frequencies ~ v cr and at early times < t x . 
However, the detailed modelling on the cyclotron emission is not important, because the flux 
is suppressed and determined by the self absorption. (We will discuss the self-absorption at 
the end of this subsection.) On the other hand, the random Lorentz factor r ) VR of the electrons 
corresponding to the optical frequency u R is relativistic. Since the electron distribution N{j e ) 
around j VR (and above it) is determined by the distribution of injected electrons at the shock, 
which are relativistic, and by the synchrotron radiation cooling, we apply the conventional 
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synchrotron model to estimates the light curve of optical flashes. 

The observer time t = (1 + z)R/ey 2 is proportional to R, because the Lorentz factor of 
the shocked shell during the shock crossing 7 X ~ mm[r],r] c ] is constant. By using the shock 
jump conditions, one finds that the scalings before the crossing time t x — (1 + z)R x jc r )\ 
are given by oc t~ 2 , oc t~ 2 and N e oc t in the thick shell case, and oc t -2 , p 3 oc t~ 3 
and N e oc t 1 / 2 in the thin shell case. The scalings of the spectral characteristics at t < t x are 
u nhr oc t _1 , z/ c r oc t and F vmax ^ r oc t° for thick shell and z/ mjr oc v c r oc £ and F v ^ max ^ r oc t -1 / 2 
for thin shell. The scalings of u Cyr and F v ^ max ^ r themselves are not correct, but the optical 
flux estimated by these scalings are right. We evaluated the scalings to calculate the optical 
light curve. 

The initial shell width A is given by the intrinsic duration of the GRB, A ~ (l+z) _1 cT 
(Kobayashi, Piran & Sari 1997), the shock crossing time t x can be written in the following 
form, 

«, ~ (|)~*T, (6) 

where 7] c ~ 60 C^^El!^ A^^T^ 1 ^ and 7\ = T/lOsec. We can determine the critical Lorentz 
factor r] c from the observations of the GRB and afterglow. If we detect the shock crossing 
time t x (the reverse shock peak time), the Lorentz factor during the shock crossing time r y x 
can be estimated from eq. (6). 

The spectral characteristics of the reverse shock emission at t x are related to those of 
the forward shock emission by the following simple formulae 1 (Kobayashi &Zhang 2003), 

V m ,r(t x ) ~ —V m ,f(t x ), V c ,r(tx) ~ V c j(t x ), F v>rnaXjr (t x ) ~ F^max, f(t x ). (7) 

Using eqs. (1), (2) and (3), we get 

VmAtx) ~ 5 -° x 1014 C 1/2 4,-i^ 2 E52 /2 ^vbZT Hz (8) 
u c , r (t x ) ~ 1.5 x 10 12 C 3/2 e- B %lE^A; 2 CR Z (9) 
V,^x) ~ 3.0 d^Ce^E^Al'^H-], Jy (10) 



1 In this paper, we assume that es, e e and p are the same for two shocked regions. Some recent works 
(Zhang, Kobayashi & Meszaros 2003; Kumar & Panaitescu 2003; Coburn & Boggs 2003) show that es,r 
might be larger than cbj where the subscripts Y and T indicate reverse and forward shocks, respectively. 
In such a case, the formulae are replaced by eqs (3)-(5) in Zhang et al. (2003). Though the power F v ^ max ^ r 
increases by a factor of (cb^A-H,/) 1 ^ 2 compared to a case of ts.r — e B.f, the cooling frequency v c decreases 
by a factor of (e Bt , r /e B j) 3//2 - This results in a dimmer optical reverse emission at the peak time by factor of 
( e s,r/ e B./) -1 ^ 4 - Note that in the ISM model, the reverse shock emission is in slow cooling regime and that 
£ s,r > gives a brighter reverse shock emission (see Zhang et al. 2003). 
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where r] 2 = i]/100 and t x ,* = t x /50sec. The reverse shock emission is in the fast cooling 
regime v cr < v mr during the shock crossing. So for t <t x , the optical flux from the reverse 
shock increases as oc t 1 ! 2 for thick shell case (vr < v m ,r(tx)) (Chevalier & Li 2000), and 
oc t( p_1 )/ 2 for thin shell case {ur > Vm >r (t x )). Since the Lorentz factor of the forward-shocked 
material is also constant during shock crossing, the optical emission from the forward shock 
evolves as t 1 ^ 2 at t < t x . But this component is usually masked by the reverse shock emission. 

At the shock crossing time t x , the optical flux reaches the peak F VRi1 . ~ F v ^ max>r 

(**Av)- 1/2 - 

F VR ,{t x ) ~ 160 d~\^e- B l\Ft A^v-rI'% %f mJy (11) 

Generally, using the time dependences of the spectral characteristics as well as eq. (7), 
we can relate this peak flux with the optical flux of the forward shock at the break time, 

/„2 \ 2-o , \ (2-3o)/4 

F, R Atx)-{^-) (j^) F URj {t mJ ). (12) 

where a — p if v m ^{t x ) is below the optical band, and a = 1 if it is above. Detection of 
the peak (t x , F VRj1 .(t x )) and the break (t m j, F UR j(t m j)) will give a constraint on the initial 
Lorentz factor rj. A similar recipe has been proposed by Zhang et al. (2003) for the ISM 
case. 

Synchrotron self-absorption would reduce our estimate (11) of the optical flash if it is 
optically thick. A simple way to account for this effect is to estimate the maximal flux 
emitted by the shocked shell material as a blackbody (Sari & Piran 1999b; Chevalier & 
Li 2000). The blackbody flux at the optical band vr and at the peak time t x is given by 
F Vr ,bb ~ 16 ( 5 ^ 4 d~ 2 v R ^ e^^A~ x E%ft t 7 J* Jy. The self-absorption frequency at the peak 

time is v sa ~ 1.0 x 10 14 Al^E^^ 1/ ' 3 i x f/ 6 H z - Considering that the peak time t x is larger 
than the burst duration T, the synchrotron self absorption will not affect the peak flux of 
an optical flash significantly for long bursts with T > 8 u^ 5 A^ 5 El^r]^ 5 sec. However, 
the self-absorption is important at early times (< t x ) in the wind model. The optical light 
curve should initially behave as t 5//2 right after the GRB trigger because of the absorption 
and it turns into ~ t 1 ^ 2 later (see Fig 1). The break time is t ~ 5 v^J 2 Alj^Elj^rf^^t^^ 
sec. 



4.3. Angular Time Delay Effect 



The cooling frequency u c ^(t x ) is well below the optical band for our typical parameters, 
the optical emission from the reverse shock should "vanish" after the peak. However, the 
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angular time delay effect prevents abrupt disappearance. The optical light curve at t > t x 
is determined by off-axis emissions. In the local frame, the spectral power is described by 
a broken power law with a low and high frequency indices —1/2 and — p/2 and the break 
frequency ~ ^ m ,r(^x)/7x- As we see higher latitude emissions, the blue shift effect due to 
the relativistic expansion becomes smaller. The blue shifted break frequency passes through 
the optical band vr at time ~ t x h> m ^{t x ) /vr. The optical flux initially evolves as t~ 5 / 2 and 
decays as £~(p+ 4 )/ 2 after the passage (Kumar & Panaitescu 2000b). 

Since the forward shock emission decays slower t -1 / 4 , it begins to dominate the opti- 
cal band at t trans ~ (t 4a /t m ,/) 1/(4a " 1) [i^ fl ,r(tx) / '-^,/(W)] 4/(4a-1) where we assumed that 
the optical reverse shock emission decays proportional to t~ a . For our typical parameters, 
v m , r (t x ) comes to the optical band. Assuming a = 13/4 (p=2.5) (see fig 1(a) in which a ~ 3 
is applied), we get t trans jt x ~ 2.5 C^^E 1 ^' 'A* 1 ^^ 1 ^ 3 ^^ 6 . The optical emission from the 
reverse shock drops below that from the forward shock within a time scale several times of 
the peak time t x . 



4.4. Duration 

If the fireball ejecta is collimated in a jet with an opening angle 9j, the duration of the 
reverse shock emission is t ang ~ (1 + z)9 2 -R x /c. The angle 6j might be determined from 
a jet break time tj ~ {0j^ x )H x of the forward shock emission (Rhoads 1999; Sari,Piran & 
Halpern 1999), even though a jet break in the wind model may not be as clear as that in 
the ISM model (Kumar & Panaitescu 2000a; Gou et al. 2001). 

^~(Mi> i/2 ~35 (_ii_y /2 (_i_y /2 min (13) 

At a low frequency v < z/ Cjr (t x ), the observer receives photons from the fluid element 
on the line of the sight until a break frequency u cut , which is equal to f Cjr (t x ) at t x , crosses 
the observational band. We now consider this time scale. After the reverse shock crosses 
the shell, the profile of the forward- shocked wind begins to approach the Blandford-McKee 
(BM) solution (Blandford & McKee 1976; Kobayashi, Piran & Sari 1999). Since the shocked 
shell is located not too far behind the forward shock, it roughly fits the BM solution. A given 
fluid element in a blast wave evolves as 7 oc t~ 3 / 8 , p oc t~ 9 / 8 and e oc t~ 3 / 2 . Assuming that 
the electron energy and the magnetic field energy remains constant fractions of the internal 
energy density of the shocked shell, the emission frequency of each electron with 7 e drops 
quickly with time according to v e oc -B77 2 oc t -15 / 8 3> 1) or v e oc B^y oc t~ 9 / 8 (7 e ~ 1). 
Using the scaling for the cyclotron emission, v cwt passes through the observational frequency 
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v at 

/ ~ 70^-4/3.-4/3 p 4/9, -16/9 ,-8/9 .13/9 . /-, ,x 

where z/ 10 = z//10GHz. Though at the radio band and early times, self absorption signifi- 
cantly reduce the flux, tdur = niax[t anfl , t cut \ can give a rough estimate of the duration of the 
radio reverse shock emission. 



5. Case Studies 

GRB 990123: The optical flash (Akerlof et al. 1999) and radio flare (Kulkarni 1999) 
associated with this burst are explained well by a reverse shock emission in the ISM model 
(Sari & Piran 1999b; Kobayashi & Sari 2000). The basic parameters of this burst include 
(e.g. Kobayashi & Sari 2000 and reference therein) E 52 ~ 140, z ~ 1.6, Ti ~ 6.3 and t x ~ 50 
sec. The wind model predicts a flatter rising of t 1 / 2 , £( p_1 )/ 2 or t 5 / 2 compared to t 3A evaluated 
from the first two ROTSE data, and a steeper decline of t -5 / 2 or /-(p+ 4 )/ 2 compared to the 
observations after the peak t~ 2 (Chevalier & Li 2000). The optical reverse shock emission is 
expected to be overtaken by that from the forward shock at t trans ~ hA* 1 ^^ 1 ^ min, but 
the observed flash decays as a single power law of t~ 2 until it falls off below the detection 
threshold at ~ 11 min. Using the jet break time tj ~ 2 day (Kulkarni et al 1999), we obtain 
tang ~ 50 min and t cut ~ 7Ae B 4/3 2 A; w/9 hr. The reverse shock emission should disappear 
well before the radio flare at ~ 1 day. We conclude that the wind model is inconsistent with 
the observations. 

GRB 021004: In the ISM model, it was shown that the major bump observed in the 
afterglow light curve around ~ 0.1 day after the burst could be explained by the passage 
of the typical frequency of the forward shock emission through the optical band, and that 
the early time optical emission is a combination of reverse and forward shock emissions 
(Kobayashi & Zhang 2003). In the wind model, the optical light curve of the forward shock 
emission initially behaves as t -1 / 4 and decays steeper as /-( 3p - 2 )/ 4 after the typical frequency 
crosses the optical band. The bump might be explained by the passage of v m j 2 . The basic 
parameters of this burst are (e.g. Kobayashi & Zhang 2002 and reference therein) E 52 ~ 5.6, 
z ~ 2.3 and Ti ~ 10. The critical Lorentz factor is i] c ~ 60A^ 4 . Assuming t m j ~ O.lday 
and F VR j(t m j) ~ lmJy, we obtain e e ~ 0.11 and e B ~ 0.085 from eqs (4) and (5). Since 
these are the typical values obtained in other afterglow observations (Panaitescu & Kumar 
2002), the wind model might also fit the observational data. 



2 After the completion of our paper, wc noticed that Li and Chevalier (2003) gave a detailed studied on 
this possibility in a recent paper. 
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6. Conclusions 

We have studied the optical emissions from reverse shocks for the thin and thick shell 
cases. The differences between this model and the ISM model are highlighted in Figure 1. In 
the ISM model, the prompt optical emission increases proportional to t 5 for the the thin shell 
case or t 1 / 2 for the thick shell case (Kobayashi 2000). In the wind model, it behaves as t 1 / 2 
for the both cases. The synchrotron self-absorption is important at early times in this model. 
The luminosity increases as t 5 ^ 2 with a steep spectral index F u oc z/ 5 / 2 at the beginning. If a 
rapid brightening with an index larger than 5/2 is caught, it could be an indication of the 
ISM-type GRB. The decay index of the emission is determined by the angular time delay 
effect in the wind model so that ~ t~ 3 , while it depends on the hydrodynamic evolution of 
the fireball ejecta in the ISM model hence ~ t~ 2 . 

When we detect the peak time of the reverse shock emission, we can estimate the Lorentz 
factor at that time from eq.(6). In the ISM model, the peak time was given by a similar 
relation t x ~ ( 7x /^)~ 8 / 3 T (Sari & Piran 1999a) where r]' c is a critical Lorentz factor and 
given by eq. (7) in Kobayashi & Zhang (2003). Additionally, if we detect the break, caused 
by the passage of v m j, in the late time (~ 1 hr after the burst) optical light curve, from 
eq. (12) we can give another constraint on the initial Lorentz factor r] (see also Zhang et al. 
2003 for the ISM case). 

This work is supported by NASA NAG5-9192 and the Pennsylvania State University 
Center for Gravitational Wave Physics, which is funded by NSF under cooperative agreement 
PHY 01-14375. 
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Fig. I. — Optical light curve: Wind Model and ISM Model. Reverse Shock Emission (solid) 
and Forward Shock Emission (dashed). 



